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摘 要 : 本 文 运用 齐 次 可 列 半 马尔 可 夫 过 程 的 向 后 方程 和 向 前 方程 ， 分 别 研究 了 GI/M/1 和 M/G/1 排 
队 系 统 队长 的 瞬时 分 布 ， 首 先 得 到 了 GI/M/1 队长 的 转移 概率 的 拉 普 拉 斯 变换 满足 的 向 后 方程 
组 ， 然 后 得 到 了 M/G/1 队长 的 转移 概率 的 拉 普 拉 斯 变换 满足 的 向 前 方程 组 ， 所 得 方程 组 的 系数 
和 矩阵 都 是 拟 下 三 角 和 矩阵 ， 都 可 以 通过 和 迭代 法 进行 求解 . 
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排队 系统 是 运筹 学 一 个 重要 分 支 ， 目 前 已 经 得 到 了 广泛 的 发 展 ， 形 成 了 一 个 独立 的 研究 方 
向 ， 近 几 年 研究 成 果 参 见 文献 [1-5] 等 ， 我 国 数学 家 越 民 义 、 徐 光辉 、 田 乃 硕 等 在 排队 论 方向 做 
了 大 量 的 研究 工作 . 

关于 GI/M/1 排 队 系 统 和 MM/G/1 排 队 系 统 的 队长 问题 ， 不 少 学 者 已 经 作 了 大 量 研究 ， 但 大 
多 是 通过 杠 入 马尔 可 夫 链 对 它们 进行 研究 ， 参 见 文 献 [6,7] 等 ， 而 本 文 是 运用 齐 次 可 列 半 马 尔 
可 夫 过 程 的 向 后 方程 和 向 前 方程 ， 分 别 对 GI/M/1 和 M/G/1 排 队 系 统 队长 的 瞬时 分 布 进行 研 
究 ， 分 别 得 到 了 GI/M/1 队长 的 转移 概率 的 拉 普 拉 斯 变换 满足 的 向 后 方程 组 ， 和 M/G/1 队长 
的 转移 概率 的 拉 普 拉 斯 变换 满足 的 向 前 方程 组 . 而且， 由 于 所 得 方程 组 的 系数 矩阵 都 是 拟 下 三 
角 和 矩阵， 所 得 方程 组 都 可 以 通过 迭代 法 进行 求解 . 


2 章 次 可 列 半 马 尔 可 夫 过 程 的 向 前 向 后 方程 


定义 1l8] SEA-WRE, X = {XX(t,w), 0 < t < 7} 称 为 齐 次 可 列 半 马尔 可 夫 过 程 ， 如 
RWE: 

1) m =0, m =inf{t > 0, X(t) # X(0)}, mai = inf{t > m, X(t) Æ X(m)},n > 1, m ÎT; 

2) P{X(Tn41) € A| Fr} = P{X(Tn41) € A| X(m)}, 1 =0,1,2,---; 

3) X(t,w) = X(m,w), Tn St < my, n=0,1,2,---. 
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WEBWi,jEeERteR,, > 
ha (t) := P(X(t) = j, t < 1 | X(0) =), P(t) := P(X(t) =5|X(0) = i), 
Q(t) := P(X (Ta) =j, m S tIX(0) =i), Qi) := P(r < t|X(0) =i), 
Gij (t) = P(Ta41 —™m < t| X (Ta) =i, X(Tny1) =j) ij 


对 应 的 Laplace 变换 记 为 
hij (A) = f e hi; (t)dt, P(A) := f e NP; (t)dt, 
0 0 


Qi) = I e*Q(dt), QWA := | e^ QM (at). 
ERI (P(t), ij € E, t e€ RR;} 是 下 述 方程 的 最 小 非 负 解 


P(t) = hyi + Y [ Qu(ds)Pej(t-s), ij EE, XER o 
k£i 


E2] {Pj(), ij € E, te Ry} 是 下 述 方程 的 最 小 非 负 解 


Pij (A) = hag (A)big + >》 Qik (A) Peg (A); i,jEE, cRy, (2) 
kżi 
方程 (1) 和 (2) 都 称 为 齐 次 可 列 半 马尔 可 夫 过 程 的 向 后 方程 . 
定理 3 外 如 果 存 在 满足 等 式 @ij( 和 ) = (WAYTARA MW Âu MP5), ij E 
E, à € Ry} 是 下 述 方程 的 最 小 非 负 解 


P(A) = hij (A) 6:3 sF NO PRAA T QARA), i, j E E, 入 E Ri, (3) 
kj 


方程 (3) 称 为 齐 次 可 列 半 马 尔 可 夫 过 程 的 向 前 方程 . 


3 ”将 齐 次 可 列 半 马 尔 可 夫 过 程 的 向 后 方程 应 用 于 GI/M/1 排队 系统 


定义 200 所 谓 GI/M/1 排 队 系统 ， 即 顾客 逐个 到 达 ， 并 且 到 达 间 隔 相 互 独立 ， 服 从 相同 
的 分 布 B(x)， 按 照 先 到 先 服务 的 原则 接受 服务 ， 只 有 一 个 服务 台 ， 每 个 顾客 的 服务 时 间 服 从 
参数 为 4 的 指数 分 布 . 

设 坟 ,表示 第 nn 个 顾客 到 达 系 统 时 看 到 系统 中 的 顾客 数 (不 包括 他 自己 )， 易 知 {7,n > 0} 是 
一 马尔 可 夫 链 . Kyo = 0 时 恰好 到 达 一 个 顾客 ，y, 为 第 n 个 顾客 到 达 的 时 间 ， 则 {Yn Yn) n > 
0} 构成 一 马尔 可 夫 更 新 过 程 . SY) = Yn, m St <n WYE) X {Yn In), 2 > 0} 伴 随 
的 齐 次 可 列 半 马 尔 可 夫 过 程 ， 而 {,n > 0} 为 {Y(t),t > 0} 的 嵌入 马尔 可 夫 链 . 

由 定义 1， 对 顾客 到 达 的 时 间 序 列 {ym} 进行 选择 ， 去 掉 拟 跳 ， 取 有 顾客 到 达 ， 过 程 的 状态 
真正 发 生 转移 的 时 刻 为 骨架 时 序列 . 在 不 引起 误解 时 ， 我 们 把 这 个 序列 也 记 为 : {yn}. 

由 半 马 尔 可 夫 核 的 定义 可 得 


Qilt ) = P(Y (y+1) = = j, m+ — In < t| Y(n) = i). 
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jti j> 0 时 ， 由 Y(t) 的 齐 次 性 可 得 
= P(Y(%m)=%m -0<tY(0) =i) 


= [ Pta 离开 了 i 二 1 一 7 个 人 }P{y € do} 
0 


f(r) ya 
Sr 
显然 要 求 i 十 1 一 7 > 0， 即 7 <i 十 1 上 式 才 有 定义 . 474i, 7 = 0 时 
i+1 i 
ri(t) = 一 : Qio(t) =Q; (t) — Ss Qi; (t) = Qi(t) > SS Gi+1~5(t) = Bit) = > gk (t), 
jH#i,j>o0 j=1,7#1 k=0, kl 
其 中 Qi 区 = P(y <tlY(0) = 引 ， 因 为 到 达 间 隔 服 从 相同 分 布 B(t)， 与 起 点 所 处 状态 无 关 ， 
所 以 Qi(t) = B(t). 于 是 GI/M/1 排 队 系 统 的 队长 Y(t) 对 应 的 半 马 尔 可 夫 核 为 
0 go(s) 0 0 
rı(s) 0 g(s) 0 
[Qut] = | 72(s) g2(s) 0 go(s) .… |, 
r3(s) q3(s) (s) 0 


上 述 结果 和 文献 [10] 中 所 得 GI/M/1 排 队 系统 对 应 的 半 马 尔 可 夫 核 的 结果 是 一 致 的 ， 只 是 去 掉 
了 拟 跳 的 情形 .Y(t) 对 应 的 半 马 尔 可 夫 核 对 应 的 矩阵 是 拟 下 三 角 的 ， 用 向 后 方程 容易 求解 . 

S Palt) = P(Y(t) =j|Y(0) = 让， 可 得 下 述 定理 . 

定理 4 对 于 GI/M/1 排 队 系 统 ，{Pj(t), ij E€ E, te R41} 是 下 述 方 程 的 最 小 非 负 解 


Po;(t) = (1 — B(t)) ðo; +/ Pi;(t — s)B(ds), 


2 二 1 
P(t) =e**(1- BO)6s+ 们 [e U ehs Pi, 1m (t — 8) B(ds). 
m=0,mA1 


证 明 对 于 G7/M/1 系 统 ， 可 得 
ho(t) := hoo(t) = P(Y (t) =0, t < um |Y (0) =0) =1- Bit). (4) 
Si AON, A PMT waz 


hi(t) := hult) = P(Y (8) = j, t < mlY(0) = i) =e"*(1- BY), (5) 
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f Qio(ds)Po;(t = 5) 一 | ri (ds) Po; (t = s) 
0 0 


i 


= f mit-aad9- S f mld- 








k=0,k#1 
‘| S (ust L (us) 
= J [ 2 k! e Hs _ 5y LT emre] P(t ~ 8) B(ds) 
0 k=0,k41 "7 Er 
t (ps)itl us 
二 f G@+D!° H Po;(t — s)B(ds). 6) 
对 任意 的 i， 下 式 成 立 
i+] Z a 
> i Quelds)Pa(t — 8) = > J qi+1-k(ds)Pk;(t — 8) 
EAS k=1,kži 
i t 
— Ss; f Qm(ds)Pi+1—m,j(t — 8) 
m=0,mA71 
a t (ps)m 
i, m! er Piti-m,j(t — s)B(ds). (7) 
由 定理 1 和 (4)-(7) 式 ， 可 得 本 定理 成 立 . 


4 
OO Oo 
n= f nla), =L A= f ad, ieN, i#1, 
0 0 


定理 5 对 于 GI/M/1 排 队 系 统 ，{Pi;( 和 ), ij E E, te Ry} 是 下 述 方程 的 最 小 非 负 解 


Poi(A) = ho(A)6oj + go(N)P1;(N), 

Pi;(N) = ja(A)617 + ri (A)Pog(A) + go(N) P23 (A), 

Poj(A) = h2(N)62; + r2(A)Poz (A) + a2 (A) Pig (A) + go (A) Ps;(N), 

Psi (A) = Pa(A)63j + ra (A) Poz (A) + 93(A) Pig (A) + g2 (A) P25 (A) + G0(A) Pa (A), 


Sio 


PI MOSHA D mAP), 


证 明 ”由 定理 2 和 GI/M/1 半 马尔 可 夫 核 对 应 的 矩阵 可 得 
Poj(A) = ho(X)éo; + $ Qor (MN) Pas (A) 
kz0 


0+1 
ho(A)603 + XO do+1—K(A)Peg(A) = po(A)5oj + G0(A) Pr (A). (8) 
k#0,k=1 
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P(A) = 访 (A)65 + 》， Qik(A)Pej(A) 
k#i 
i+1 
hi(A)biz + Qio A) Poz (A) + > Qik (A) Pij(A) 
k#i,k=1 


i+1 
hi(A)big + 7i(A)Pg(A) + >, giga—e(A) Pag (A) 
kfi,k=1 


hi(A) big + 7e(A)Pos(A) + YO Qn(A)Pit1—m,j(A)- (9) 


m#1,m=0 


i 


由 定理 2 和 (8), (9) 式 ， 可 得 本 定理 成 立 . 
由 定理 5 可 知 ，{Pi;( 和 ), ij € E, t € BR+} 满 足 的 向 后 方程 组 对 应 的 系数 和 矩阵 是 拟 下 三 角 
的 ， 可 以 用 迭代 法 进行 求解 . 


4 ”将 齐 次 可 列 半 马尔 可 夫 过 程 的 向 前 方程 应 用 于 M/G/1 排队 系统 


定义 309 ”所谓 M/G/1 排 队 系 统 ， 即 顾客 逐个 到 达 ， 并 且 到 达 的 间隔 服从 参数 为 和 的 指数 
分 布 ， 按 照 先 到 先 服务 的 原则 接受 服务 ， 只 有 一个 服务 台 ， 每 个 顾客 的 服务 时 间 独 立 同 分 布 ， 
服从 分 布 4(z). 

设 X 表示 第 n 个 顾客 离开 系统 时 系统 剩余 中 的 顾客 数 (不 包括 离开 的 这 个 人 )， 易 知 {Xn 
> 0} 是 一 马尔 可 夫 链 . Bro = 0 时 恰 有 一 个 顾客 离开 ,7 为 第 n 个 顾客 离开 的 时 间 ， 
则 {(X,74),n > 0} 构成 一 马尔 可 夫 更 新 过 程 . AX) = Xn, mm <t< trary MXA (Xn, 
Tahin > 0} 的 伴随 的 齐 次 可 列 半 马 尔 可 夫 过 程 ， 而 {Xn > 0} 为 {X(b,t > 0} 的 嵌入 马尔 可 
K. 

由 定义 1， 对 顾客 离开 的 时 间 序列 (7, 进行 选择 ， 去 掉 拟 跳 ， 取 有 顾客 离开 ， 并 且 过 程 的 
状态 真正 发 生 转 移 的 时 刻 为 骨架 时 序列 ， 在 不 引起 误解 时 ， 我 们 把 这 个 序列 也 记 为 ; {rn} 

由 半 马 尔 可 夫 核 的 定义 可 得 


Out) = P(X Gas) = J, n41 — Tn S t| X (Ta) =å). 
当 7 关 i i > 0 时 ， 由 X(t) 的 齐 次 性 可 得 
Qj-it1(t) =: Qi(t) = P(X(tn41) = j, mii ~ Tn < t| X(t) =å) 


= P(X(n) =j, 1-0 < t|X(0) =i) 


(ALEE ig 
= f 0 Md A(z), 
显然 要 求 7 一 i 十 1 > 0， 即 i <j 了 十 1 上 式 才 有 定义 . 

当 j 关 i, i =O 


pi(t) =: Qoz (t) = P(X(n)=j mn — 0 < t| X(0) = 0) = [ ue- aean 
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于 是 MVG/1 排 队 系统 对 应 的 半 马 尔 可 夫 核 为 


0 pi(s) pa(s) pa(s) = 
go(s) 0 9%(s) ga(s) 

[Qn] =| 9 ws) 0 gals) …|， 

0 0 qs) 0 


上 述 结果 和 文献 [10] 中 所 得 M/GV1 排 队 系统 对 应 的 半 马 尔 可 夫 核 的 结果 是 一 致 的 ， 只 是 去 掉 
了 拟 跳 的 情形 ，MV/G/L 的 半 马 尔 可 夫 核对 应 的 矩阵 是 拟 上 三 角 的 ， 这 种 情形 下 ， 用 向 前 方程 
比 向 后 方程 更 容易 求解 . 

令 


P(t) = P(X) =j|X(0) =i), P(A) = f o-Ps (t)dt, 

i( 和 A) = * ety, (dt), (=1,2,---, g(A)= * eg (dt), EN, i£l. 

pil) f pi(dt), i qa) [e a(dt), i€N, iF 
定理 6 对 于 M/G/1 排 队 系统 ，{Pj(A), ij EE, te Ri} 是 下 述 方程 的 最 小 非 负 解 


Pyo(A) = ho()5i0 + Pin (A) (ha (A) oo(Aho(A)， 

Pix (A) = hi (Abia + Pio(A) (ho(A)) pi (A) (A) 十 P(A) (h2(N)) go(A)ja(A)， 

Pa(A) = hz(N)6i2 + Pio(A) (ho (A))7"po(A)ha(A) + P(A) (ha (A) g2(A)ha(A) 
+ Pia (A) (ha (A)) T go(A)h2(A), 

P(A) = hi (A) big + Piol A) (hA T PAA) 


+ Pan-n O) (hym) mO, 


证 明 ”由 定理 3 和 M/G/1 半 马尔 可 夫 核 对 应 的 矩阵 可 得 


Po(A) = ho(N)6io + YO P(A) (he(N)) Quo(A)ho(A) 


k#¥0 


0 十 1 


hoNiiot+ YO Ph(A)(Ak(A)) oo-eHi(A)ho(A) 
k#0,k=1 


= ho(A)di0 + Pa (A) (h (A)) “go(A)ho (A). (10) 
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当 了 A ORT 
PaA) = 万 (65+》 P(N (hi(N)) Qn AA) 
kAj 
j+1 
= hj(A)6iz + Pio(A)(ho(A))” RAA 》 PAAA) "Ques (A)RG(A) 
kA~j,k=1 
了 十 
= hj(A)6i + Pio(A) (hA) PANA YO PRAA) oO 
kAj,k=1 
= h;(A)éy + Pio(A)(ho(A)) T B A)h (A) 
j 
+ JO PstimO) (hrm(N)) am (A)hi(A). (11) 


mA1,m=0 


于 是 ， 由 定理 3 和 (10), (11) 式 ， 可 得 本 定理 成 立 . 

由 定理 6 可 知 ，{Pij;( 和 ), i,j € E, t € RH 满足 的 向 前 方程 组 对 应 的 系数 矩阵 是 拟 下 三 角 
的 ， 可 以 迭 用 代 法 进行 求解 ， 而 它 满足 的 向 后 方程 组 对 应 的 系数 矩阵 是 拟 上 三 角 的 ， 不 利于 求 
解 ， 这 体现 了 向 前 方程 组 的 应 用 价值 . 
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The Application of Semi-Markov Processes in GI/M/1 and M/G/1 
Queuing Systems 


DONG Hai-ling’, HOU Zhen-ting?, JIANG Guo-chao® 


(1- School of Mathematics and Computer Science, Shenzhen University, Shenzhen 518060; 
2- School of Mathematics, Central South University, Changsha 410075; 
3- Department of Economics and Management, HIT Shenzhen Graduate School, Shenzhen 518055) 


Abstract: In this paper, the backward equations and forward equations of homogeneous denumerable 
semi-Markov processes are used to study the instantaneous distribution of queue lengths of GI /M/1 
and M/G/1 queuing systems, respectively. For the GI/M/1 queuing system, the backward equations 
of the Laplace transform of the transition probability of its queue length are obtained. For the M /GI/\ 
queuing system, the forward equations of the Laplace transform of the transition probability of its queue 
length are obtained. These backward equations and forward equations, whose coefficient matrices are 
quasi-triangular matrices, can be easily solved by using the iterative method. 

Keywords: homogeneous denumerable semi-Markov processes; GI/M/1 queuing system; M/G/1 
queuing system; backward equations; forward equations 
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